Active Absorbing State Phase Transition Beyond Directed Percolation : A Class of 

Exactly Solvable Models 
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We introduce and solve a model of hardcore particles on a one dimensional periodic lattice which 
undergoes an active-absorbing state phase transition at finite density. In this model an occupied site 
is defined to be active if its left neighbour is occupied and the right neighbour is vacant. Particles 
from such active sites hop stochastically to their right. We show that, both the density of active 
sites and the survival probability vanish as the particle density is decreased below half. The critical 
exponents and spatial correlations of the model are calculated exactly using the matrix product 
ansatz. Exact analytical study of several variations of the model reveals that these non-equilibrium 
phase transitions belong to a new universality class different from the generic active-absorbing-state 
phase transition, namely directed percolation. 
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One of the most studied models of non-equilibrium 
phase transition is directed percolation(DP)|lj] defined on 
a d-dimensional lattice, where an infected site stochas- 
tically infects its neighbours in one particular direc- 
tion. Depending on the infection probability p, the 
infection may eventually survive (when p > pc) or 
decay into the absorbing state where no site is in- 
fected. Non-equilibrium phase transitions in several 
other systems, such as reaction diffusion systemsQ, de- 
pining transitions [s^l , damage spreading synchroniza- 
tion transition 5], sand-pile modelsQ, and certain prob- 
abilistic cellular automata^, are known to be in the uni- 
versality class of DP. It has been conjectured [1| that an 
"active-absorbing phase transition governed by a fluc- 
tuating scalar order parameter" generically belongs to 
the universality class of DP. There are certain excep- 
tions, though. Particle-hole symmetry^, conservation 
of parity I IQl . and symmetry between different absorb- 
ing states[ll| lead to different universalities. Again in 
sandpile models coupling of the order parameter 
to the conserving height fields [l^ results in different 
critical behaviour. Also, conserved lattice gas (CLG) 
models [3, [3 where the activity field is coupled to the 
conserved density show critical behaviour different from 
DP. This absorbing state phase transition in the presence 
of conserved field is not well understood and most studies 
in this direction are numerical. 

In this paper, we provide an exact analytical solution 
for a model of hardcore particles on a one dimensional 
ring which undergoes an active-absorbing phase transi- 
tion as the density of particles is changed. The model is 
defined with a dynamics where a particle from an occu- 
pied site hops to the right neighbouring site if the left one 
is occupied. This restricted asymmetric exclusion process 
(RASEP) leads to a transition from an active phase to an 



absorbing state as the density of the system falls below 
i. The critical exponents of the system at the transi- 
tion point and spatial correlations have been calculated 
exactly using the matrix product ansatz (MPA) Il7| . 
Some variations of the model, where particles may hop 
to both directions stochastically, or hop to the right (left) 
only when it is followed by or more particles from 
left (right), could also be solved exactly. These models, 
which have same exponents at the transition point, form 
a new universality class of active-absorbing phase transi- 
tion different from the generic universality class, namely, 
DP. 

The model is defined on a one dimensional lattice la- 
beled by sites i — 1,2 ... L which are either vacant or 
occupied with at most one particle; corresponding site 
variables are taken Si — 0,1. A periodic boundary con- 
dition is imposed so that s^+l — Si- The dynamics of 
the system can be described as follows. A particle from 
a randomly chosen site i is transferred to the right only 
if Si+i = and Si_i — 1. This particle conserving dy- 
namics is thus equivalent to a reaction diffusion system 



110^ 101. 
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We define the activity field at site i as = Si^iSi{l — 
Si_|_i) which takes values 1 or depending on whether 
the site i is active, i.e., Si — 1 — Si-i andsi+i = 0. The 
density of active sites 



(Sl_lSj(l 



)) 



(2) 



* E-mail address: lurna.basu@saha.ac.inl 
t E-mail address: |pk.mohanty@saha.ac.in"] 



is denoted by pa in the thermodynamic limit. A config- 
uration is said to be active if there is at least one active 
site, otherwise it is called absorbing. For a system of 
N — Si particles density is p = ^. Clearly, there is 
only one configuration at p = (and at p = 1) which is 
absorbing. First let us consider the regime p < ^ where 
there are both active and absorbing configurations. To- 
tal number of absorbing configurations in this regime is 
L-N ^N^^ ^iid the rest are active. In this regime the sys- 
tem is arrested in one of these absorbing configurations 
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equation in steady state, 



FIG. 1: The order parameter pa for RASEP is nonzero for 
p > 1/2. Inset shows the decay of correlation function Tj 
defined in Eq. ((15} . 



in steady state resulting in pa = 0. For p > i, however, 
there is no absorbing configuration which corresponds to 
an active phase with fluctuating density of active sites. 
Thus, Pa can be taken as the order parameter of this ac- 
tive absorbing state phase transition occurring a,t pc — \. 



One can describe the dynamics of the model alter- 
natively in terms of the bond variables r^, which con- 
nects the sites Si and Si+i. Correspondingly, we choose 
Ti = 2si + Si+i for four possible combinations {si,Si+i). 
Note that every configuration {si\ can be uniquely trans- 
lated to {ti] and vice verse. The dynamics can be de- 
scribed in terms of Ti as a reaction system 
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The fact that Ti and r^+i have a common site Si+i puts 
certain restrictions on the allowed configurations. How- 
ever one need not bother about those as any configuration 
translated from {s^} automatically satisfies these restric- 
tions and the dynamics ^ respects the same. 

To get a steady state distribution for the reaction sys- 
tem, either Eq.((T]) or Eq.Q which has three site dy- 
namics, we generalize the formulation of matrix prod- 
uct ansatz [l^ [I3| which is commonly used for a two 
site dynamics. This generalization is different from what 
has been discussed earlier 18]. Let us describe the for- 
mulation in generic terms before using it in this specific 
problem. In MPA, first a configuration {ci, cr2, . . . ctl} is 
translated to a product of matrices by replacing each ai 
by a matrix and an ansatz is made that for a peri- 
odic system, the unnormalized weight in the steady state 
is given by 



/(cri,cr2, . . .CTl) = Tr[A„^,Aa 



(4) 



This ansatz could provide an exact solution for any 
three site dynamics if one can find matrices such that 
steady state weight Q satisfies the corresponding master 



= 



'^W{aiai+iai+2 CT-cr-_,_ia--+2)/(---CTjcrj+icri+2- 



Here Ws are the transition rates for the three site dy- 
namics. Right hand side of the master equation can be 
arranged to vanish for any generic three site dynamics if 



W{<7i(Ji+iai+2 cr-CT-+iCr^_|_2)/(...o-iCri+iCri+2---) 
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where Aa- are auxiliary matrices. Equation ([5]) is a gener- 
alization of the cancellation procedure introduced earlier 
[l^ for two site dynamics. Such a cancellation is feasible 
only when one can find matrices and auxiliaries which 
satisfy Eq.© for a specific dynamics. 

Now let us try to apply this generic scheme to the 
dynamics ^ and ([3]). In the first case ([1]), by replacing 
Si with a matrix Ag. we find that the cancellation would 
occur only if 

AiAiAo = '^AiAiAo + AiAiAo 

= AiAoAi- AiAoAi. (6) 

Note, that these algebraic relations can not be satisfied 
by non-zero scalars Aq and Ai, but there are solutions 
where An, A^ and the auxiliaries are finite dimensional 
matrices [19l|. 

Next, for the dynamics ^ with bond variables we re- 
place Ti by matrices and demand that the generic 
cancellation scheme (O should hold for this dynamics 
([3|). Then, the matrices {A",-} and auxiliaries {^t} must 
satisfy 

X:jX2Xi — A^2^i;^3 ~ X2XiX-i^ 

= —X^X2Xi + A'3Ar2A"i (7) 

A'3A^2^0 = X2X1X2 — X2XiX2^ 

— —X3X2X0 + A'3A'2A'o (8) 

It is not difficult to see that Eqs. ([7]) and ^ have 
a scalar solution Xq — 0,Xi = X2 = X3 = 1 with 
auxiharies Xq — 0,Xi = 2 = X2, X3 = |. Usually 
the solutions of MPA with one of the matrix being zero 
are not acceptable as it indicates that certain configu- 
rations are never visited in steady state. For example, 
here Xq = would mean that steady state weight is 
zero for all configurations having two or more consecu- 
tive zeros. Thus, to accept solutions with Xq = we 
must show a priori that the steady state of the above 
said configurations are in fact zero. A direct proof is 
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lengthy. Alternatively, one can prove the same using a 
mapping of the model to zero range process (ZRP) which 
is discussed later [below Eq. (PT|) ]. Since it is easy to work 
with scalars, we choose to continue with dynamics 
Let us first calculate the partition function keeping in 
mind that (i) and r^+i have a common site Si+i, and 
(ii) Si — N. The first restriction can be taken care 
of by defining ~ {si\Y\si+i). Using the above scalar 

^ ^ Then the partition function is 



solution, Y 



1 1 



Y,'l[{s,\Y\s,+,) 



(9) 



where ' reminds that the sum is restricted by Si — N. 
To evaluate the restricted sum we go over to grand canon- 
ical system which is an ensemble of L site rings, each 
having a weight where N, the number of particles in 
the ring, takes all possible values. The grand partition 
function 



where T ■ 
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Tr{T^) =\% + \^_ (10) 







and A± = f (z ± \/4 + z^) are the 
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eigenvalues of T. Average density of particles is then 



p{z) 



lim 

L — >oo 
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= A+/(A+ + A_), 
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where L oo limit has been used in the last step. A 
system with fixed density p = N / L would correspond to 
the choice of z which is consistent with Eq. (fTT|) . i.e.. 



2p- 1 

v/p(i-p)' 



(12) 



Now, let us calculate some of the observables. First, the 
order parameter. 
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(i|r|i)(i|T|o)(o|r^-^|i) 



-l-A^(A_/A+)^-2 



Pa 



where 

Pa = (2p-l)(l-p)/p 



(13) 



is the order parameter of the system in the thermody- 
namic limit L oo. As pa — ed p — Pc — 1/2, {(f>i) 
is independent of L at the critical point. However for 
p > Pc {4>i) has a finite size correction and it converges 
to Pa exponentially with L. 

Any observable can be calculated from the generic (n+ 
l)-point correlation function, 



= lim — — r-^ 



l|r|l)"(l|T 



ID 



p 



2p-l 
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For example, (</),;) = {siSi+i{l - 5^+2)) = Ci - C2 = 
(2p — 1)(1 — p)/p, which is same as Ea. ((T^ . Now, as- 
suming translational invariance, correlation between two 
active sites separated by j lattice sites, Tj — — 
can be calculated as follows. We have. 



— (l|T|l)2(l|T|0)2(0|r^-2|i)(0|r^-J-2|i) 



Pa 



1 - 



resulting in 



T,ip) = -{2p-l)' 



P-1 
P 



(15) 



Note that Tj (p) oscillates with j as shown in the inset of 
Fig. 1. Such an oscillation is expected as 0i(/)i±i = = 
0i'/'i±2 for every i. 

Let us calculate the critical exponents of this phase 
transition. Formally the correlation function is written 
as r(j) - exp{-j/l)j'^-^+'^-i\ Thus, Eq.dUD implies 
that 7y = 1 and that the correlation length ^ = (In j--^)^^ 
diverges as ^ ^ (p — Pc)^'' with v — 1. From p^ . pa is 
linear in {p— Pc) near the critical point. So, the order pa- 
rameter exponent [3—1. Again, the survival probability 
P that a single active site survives in f — > 00 limit, van- 
ishes as P = (p— pc)^ . In RASEP, the activity certainly 
survives for any density p > pc\ thus, /3' = 0. 



TABLE I: Critical exponents of DP, CDF, and RASEP 





P 




P' 


Dpyj 


0.276486 


1.504144 


1.0968 0.276486 


CDP[1J 





1 


1 1 


RASEP 


1 


1 


1 



For comparison we have listed all these critical expo- 
nents of the model along with those for other known uni- 
versality classes of active-absorbing phase transitions in 
Table L The most well-known and generic universality 
class of active-absorbing phase transition having a fluc- 
tuating scalar order parameter is DP[8|. Models where 
order parameters obey special conservation laws could 
differ from DP. One such example is compact directed 
percolation (CDP)'^ where the activity field satisfies the 
particle hole symmetry. In RASEP, the order parameter 
which is scalar and fluctuating, does not satisfy any spe- 
cial conservation law. That, it shows an active-absorbing 
phase transition different from DP is surprising. Cou- 
pling of this fluctuating order parameter to a conserved 
field, namely density, could be a possible cause. In fact, 
it is well known [Ijl that in sand-pile models of self or- 
ganized critically, the activity field (which is scalar and 
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fluctuating) is coupled to the conserved height field re- 
sulting in universality classes different from DP. 

To know if this new universality class is stable against 
perturbations we have studied several variations of the 
model by introducing stochasticity both in the direction 
and rate of particle transfer, 

110^101 and 011^101 (16) 

Naturally, here a site is called active when Si — 1 and 
either of Si±i — 0. Note that CLG in one dimension is a 
special case of Eq. (IT6)) with p — 1 where both forward 
and backward hopping of particles are allowed. Since 
this symmetric dynamics satisfies detailed balance [Tsf. all 
the allowed configurations have the same weight in the 
steady state . For generic p 1, however, the dynam- 
ics (|16l) does not satisfy detailed balance. To obtain the 
exact steady state distribution for arbitrary < p < I 
we use MPA for three site dynamics described in this 
paper. Here again, configurations with two or more con- 
secutive zeros are never visited in the steady state, and 
all other configurations are equally probable. The tran- 
sition occurs at the critical density pc — ^. The density 
of active sites in the active phase (p > pc) is given by 
Pa ~ 2(1 — p)(2/9— 1)//9, which vanishes linearly as p ^ Pc 
resulting in /3 = 1. Other critical exponents v = rj = 1 
are found to be the same as that of RASEP. 

It is worth mentioning that a transformation 1 ^ of 
(fT6|) which leads to a dynamics 

001^010 and 100^010 (17) 

also shows a transition at pc = 1/2, with order parame- 
ter Pa = 2p(l - 2p)/(l - p) for p < 1/2. The critical be- 
haviour here, as expected, is same as that of RASEP. An 
interesting variation is when both the dynamics p6|l and 
(|17p are present. In this case we have only two absorbing 
states {101010 . . .} and {010101 . . .} which are symmet- 
ric. This may lead to different critical behaviour [TJ], as 
supported by the numerical studies of these models for 
P = 1 20]. 

In another variation, a particle from an occupied site 
hops to its right only when it is followed by p occupied 
sites from its left. Thus, Eq.([T]) is a special case with 
p = 1. For finite /i = 2, 3 . . . the dynamics are 

p = 2 : 1110 1101 

p^i : 11110 -> 11101 ... etc. (18) 

To use MPA for these dynamics, we have extended our 
formulation for {p 4- 2)-site dynamics. The exact results 
there show that this class of models with p > 2 undergo 
an active absorbing phase transition which belong to the 
same universality class as the system with p = 1. Details 
of these calculations will be published elsewhere. Here we 
show a mapping of these models to the ZRP 21] which 
simplify calculations of some of the observables, such as 
Pa and its distribution. 

The ZRP is defined on a periodic one dimensional lat- 
tice with the following dynamics; a single particle from 



a randomly chosen site (or box) hops to one of its neigh- 
bour, say the right one, with rate u{n) which depends 
on the number of particles n in departure box. To map 
Eq. (fT8|) to ZRP, we define the vacant sites as boxes, and 
the number of uninterrupted sequence of Is to the left 
of a vacant site as the number of particles in that box. 
Thus there are N particles which are distributed among 
M — L — N boxes. Now, dynamics HH]) just transfers a 
particle from a box to its right if the departure box has 
more than p particles. Thus, 

u{n)^e{n~ p), (19) 

where 9{x) is the Heaviside theta function. Clearly, when 
the particles per box g — jj < p^ the system has at 
least one configuration where every box contains < p 
particles. Such configurations are absorbing and the 
system is arrested in one of them in the steady state. 
Thus the critical density is Qc — p, which corresponds to 

Notice that the steady state weight in the active phase 
p > Pc has product measure : 

P{ni,n2 . . . hm) ^ f{ni)f{n2) . . . /(^m), (20) 

where Ui is the number of particles in box i and function / 
is to be determined such that Eq. ((20| satisfies the master 
equation in steady statef2i]. In this case, the rate of 
transfer is independent of number of particles resulting 
in 

f{n)^e{n + l-p). (21) 

Thus, in steady state, all configurations with every box 
containing p or more particles are visited with equal 
probability and all other configurations, which have at 
least one box containing less than p particles, are never 
visited. In particular for RASEP ([1]), configurations hav- 
ing two or more consecutive zeros are not allowed in 
steady state (as claimed earlier). 

Partition function of the system, in this case, is just 
the total number of configurations where N particles are 
distributed in i — boxes such that each contains at 
least p particles. 

Note that every configuration {s^} has L translationally 
equivalent configurations, whereas in ZRP it has only 
only L ~ N equivalent ones. This raises a multiplicative 
factor L/{L — N) to the steady state weight of every con- 
figuration. We have ignored this factor in Eg. iP^ and in 
further calculations as it does not affect the observables. 
This mapping allows the calculation of fluctuations in 

the number of active sites Na — "Pi — Ef=i ^("-j ~ 

p). In a system of N particles distributed among M 
boxes probability of finding Na boxes which have more 
than p particles with a restriction that every box con- 
tains at least p particles is given by a hypergeometric 
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distribution, 

P{Na) 



1 ^N-^i{L-N)-l ^L-N 



(23) 



Mean and variance of this distribution are related to the 
order parameter and its fluctuation respectively as 



Pa 



= lim 



L p - {p^ - - p) 



1 



hm -{{Nl) - (N^y) 

L— i-oo L/ 

Pa' 

p-{p^l){l-p) 



(24) 



(25) 



Equation ([M]) provides an exact expression of pa for the 
generic model, which vanishes linearly as p approaches 
Pc = ifji- From Eq. p5|) it is clear that the fluctua- 
tion vanishes quadratically as p —^ pc- Contrary to other 
known continuous transitions, here the transitions are 
not associated with diverging fluctuation of the order pa- 
rameter but it is associated with a diverging correlation 
length. 

Of course, one can extend these models to incorporate 
particle transfer to both directions. A model of ZRP hav- 
ing a generic threshold /i and unbiased particle transfer 
has been studied earlierfil] in d-dimensions. Character- 
istic critical exponets, in the context of sand-pile models, 
have been discussed. 



In conclusion, we have introduced a class of models 
in one dimension where a particle can move to a vacant 
neighbouring site in one direction only if it is followed 
by p number of particles in the other direction. We ex- 
tended the matrix product ansatz to generic three site 
dynamics and apply the formalism to the simplest ver- 
sion of the model with p = 1 to get the exact steady state 
distribution. We show that these models undergo a con- 
tinuous active-absorbing phase transition when density 
of particles is decreased below pc — p/{l + p). Inter- 
estingly the fluctuation of the order parameter here does 
not show any divergence at the transition point, whereas 
active sites are found to be correlated within a length 
scale ^ which diverges as critical density is approached 
from above. Critical exponents of this active-absorbing 
phase transition, which are calculated analytically, are 
found to be different from the generic universality class, 
namely directed percolation. We argue that the fluctuat- 
ing scalar order-parameters in these models are coupled 
to the density fleld which is conserved, which could be 
a possible reason why these models differ from the DP 
universality class. 
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